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ABSTRACT 


A series of tests were conducted at JPL to measure the mass flux in 
the far field of a nozzle plume in a high vacuum with emphasis on the back 
flow region. Existing theories to predict the far field of a plume are not 
adequate for large angular departures from the plume axis. The measure- 
ments presented in this report provide fairly accurate data for off-axis 
angles as large as 140° (i.e. , in the back flow region). This region, since 
it is well behind the exit plane, is of particular interest to those concerned 
with instrument contamination. Usually sensitive spacecraft surfaces are 
located in the region affected by the back flow. 

The tests, which utilized five different nozzles, were performed at the 
JPL Molsink facility. Parameters such as expansion ratio, throat diameter, 
nozzle lip shape, and plenum (chamber) pressure were varied. Carbon 
dioxide and nitrogen gases were flowed and mass flux measurements were 
taken using quartz crystal microbalances in as many as nine different loca- 
tions relative to the test nozzle. 

The tests have resulted in a large matrix of data that were correlated 
and compared to the Hill and Draper flow prediction theory. These tests are 
a continuation of earlier attempts to provide quantitative data, the results of 
which were previously published in two JPL reports. 

Several conclusions with respect to the effect of nozzle and gas param- 
eters on the amount of back flow mass flux are offered, and it was demon- 
strated that gaseous mass fluxes, which are not predictable by present 
theories, are encountered in the region behind the nozzle exit plane. This 
knowledge is significant if materials incompatible with the gaseous exhaust 
products are used in this region. 
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I. INTRODUCTION 


Plumes from small attitude propulsion engines as well as from 
spacecraft main engines have in the past been suspected of impinging on 
spacecraft surfaces and of causing anomalous behavior of the instrumenta- 
tion. Existing plume prediction methods fail to give a satisfactory account 
of the presence of plume gases in the far upstream regions (large plume 
turning angles) where most instrumentation is located. The problem is 
increased if the Reynolds number of the internal flow becomes sufficiently 
small that the boundary layer occupies almost the entire flow field, since it 
is the gas in the boundary layer which turns beyond the predicted limits. 

This phenomena, referred to as boundary layer expansion around the 
nozzle lip, has been treated in the literature by several investigators (see, 
for example, Reference (1)). In some instances, as in the cited example, 
semi -empirical methods are proposed as an expedient to account for the 
boundary layer effects in the plume far field. However, the criterion of 
comparison proposed by these methods has been based exclusively on numer- 
ical calculations which, although lengthy, fail to include in the modeling 
the subsonic part of the boundary layer. It is, therefore, highly desirable to 
obtain accurate experimental data on which to base an empirical theory. The 
topic of this paper is to present a collection of this type of data obtained 
recently at the Jet Propulsion Laboratory. 

The series of tests reported herein are a continuation of earlier 
attempts to provide quantitative data on this subject, the results of which 
have been published in two JPL reports (References (2) and (3)). These 
tests, which utilized five different nozzles, were performed at the JPL 
Molecular Sink, high vacuum facility (Molsink). Carbon dioxide and nitrogen 
were selected as representative gases, and quartz crystal microbalances 
were used in multiple locations to measure the mass flux in the far field of 
the plume. 
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The results presented in this paper are accurate within a few percent 
and do accurately reflect the physical phenomena which occurs when gases 
from nozzles with large boundary layers expand into a vacuum. The results 
were reduced to a form amenable for comparison with existing theories, and 
it was demonstrated that significant gaseous mass fluxes are encountered in 
the back flow regions where existing theories are not applicable. These 
findings are particularly meaningful to those concerned with the effect of 
rocket exhaust on sensitive surfaces of the spacecraft or on scientific pack- 
ages that could be directly or indirectly affected by the presence of small 
traces of exhaust gases. 


II. EXPERIMENTAL 


A. F acility 

One of the major problems associated with the study of vacuum nozzle- 
plume flow fields has been the sensitivity of the plume far field to the chamber 
recirculation effects. Since the measurements that this work was concerned 
with were specifically in the far field, the Molsink facility was selected, as 
it is an ultra-high vacuum facility capable of cryopumping injected gases at 
a very high rate. This facility consists of the vacuum chamber and the asso- 
ciated cryogenic systems (see Figure 1). The vacuum chamber encloses two 
other chambers: the inner liner and the molecular trap. The inner liner is 

filled with liquid nitrogen that acts as a massive heat sink. The molecular 
trap (moltrap), the innermost chamber, is a sphere approximately 3 m (10 ft) 
in diameter, maintained at a temperature between 10 and 15°K with gaseous 
helium. The walls of the moltrap are wedge-shaped, resembling an anechoic 
chamber, with a total surface area of approximately 186 m 2 (2000 ft 2 ). The 
liquid nitrogen is supplied to the inner liner from a central tank, and vented 
to the atmosphere. The moltrap is cooled by a manifold of tubes within which 
gaseous helium at approximately 7°K circulates. The helium is kept at this 
very low temperature by a refrigerator located adjacent to the chamber. The 
behavior of a rocket plume inside the Molsink can be described by comparing 
the flow field both in space and inside the chamber (see Figure 2). The rocket 
exhaust in space expands freely in an almost radial flow. If a hypothetical 
perfect sink surface were to enclose such a plume, the flow field enclosed by 
such a surface would be identical to the one . experienced in space. This is 
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the working principle of the Molsink chamber. However, since the walls of 
the chamber are not perfect sinks, a small molecular reflection occurs. This 
reflection results in a recirculation effect that could degrade the space simu- 
lation. Since the reflection coefficient depends on the vaporization rate of 
the gases at the wall temperature and diminishes as the wall temperature is 
decreased, the vaporization rate of the gases used in these measurements, 

N 2 and CC >2 , is for all practical purposes negligible at 10°K. More informa- 
tion on the facility is available in References (4) to (6). 

B . Test Apparatus 

In order to assess the effect of nozzle parameters on the plume far field 
characteristics, five different nozzles were fabricated. Each nozzle was pro- 
vided with its own cylindrical plenum as shown in Figure 3. The plenum was 
connected to a length of heavy wall tube 2.13 m (7 ft) long as shown in the 
figure. For each of these five nozzles a thermocouple and a pressure tap 
were installed on the plenum, and a coaxial heater was introduced into the 
connecting tube to provide adequate thermal environment to the tube-plenum- 
nozzle assembly. The terminals of the plenum pressure tap, heater, and 
thermocouple were rigidly fixed to the plenum. At the other end of the tube 
the instrumentation was inserted through leak-proof feed-throughs which were 
part of the tubing assembly. In this manner, the tube-plenum-nozzle assem- 
bly and its instrumentation constituted an independent unit. 

The fluid wetted surfaces of each of these units were reproduced as 
close as possible in order to maintain the same conditions upstream of the 
nozzle throat. The units were installed into the chamber by introducing the 
entire set through the Molsink top door. The axis of each of the five nozzles 
was positioned in a pentagonal pattern as shown in Figure 4. They were held 
in place by a perforated aluminum plate and its feed-through inserts. The 
aluminum plate assumed the role of the upper door, and an additional guide 
for each of the tubes was provided by a perforated plate of micarta that fitted 
the aperture at the top of the moltrap. With this arrangement each of the 
five units could be moved up and down and still maintain a very small toler- 
ance for lateral departure. The tubes were fed through tight neoprene fittings 
which were installed in each hole in the aluminum plate. Lubrication for the 
tube motion was provided by DC -11, a low outgassing vacuum oil lubricant. 

The end of each tube outside the Molsink was sealed with a hand valve. The 
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gases were obtained from bottles located on top of the chamber and injected 
into the tube -plenum-nozzle unit by means of a combination of regulator and 
pressure gauges. The thermocouples were connected into the data system 
and the pressures were recorded from pressure transducers connected to 
the pressure tap terminals. Table I summarizes the characteristics of 
these nozzles . 

C . Instrumentation 

Quartz crystal microbalances (QCM) were used to measure the mass 
flux in the plume far field, which will be defined herein as all distances 
greater than 20 exit diameters from the nozzle. 

A QCM consists of an electronic oscillator whose resonance frequency 
is stabilized by the piezoelectric effect of a quartz crystal. The resulting 
resonance frequency depends on several parameters, but if one fixes the 
oscillator circuit constants and polarization voltages, the specific modes of 
crystal vibration will depend only on the orientation of the cut plane with 
respect to the crystal axes. Depending on the angle of the crystal cut, the 
precise resonance frequencies will be a function of both the mass deposited 
on the surface of the crystal and the temperatures. If the crystal experiences 
a AT variation in temperature and a AM mass variation, the frequency shift 
can be expressed as: 


Af = C m AM + C t AT 

where C M and C T depend on the temperature and cut angles of the crystals 
(see Nomenclature). If a cut angle is chosen such that C T = 0 for some 
range of temperatures, then Af = C^AM and the crystal can be used as a 
delicate microbalance to detect and measure small masses deposited on the 
surface. Since for a considerable change in temperature and cut angle of the 
crystal selected, the mass coefficient does not vary more than ±5%, one can 
use, for all practical purposes, the expression 

(D 

AM = 4. 30 X 10 ' 7 =y 

F 

c 

where M = mass deposits in g/cm 2 , f = frequency shift in Hz, and F c = 
resonant frequency in MHz. 
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An illustration of the crystal, along with the location of one of the 
electrodes and the thickness vibration mode, can be seen in Figure 5. Addi- 
tional information related to QCM theory can be found in Reference (4). 

The crystal as described in the above paragraph will trap a gas by 
condensation if its surface is exposed to the flow and its temperature kept 
well below the condensation temperature of the oncoming gases. Since nitro- 
gen and carbon dioxide were used for this set of experiments, temperatures 
as low as 15°K were used in order to maintain a very low evaporation rate. 
Thus, a particular crystal cut with a large coefficient and a low Crp coef- 
ficient at temperatures below 15°K had to be determined. It was found that 
an AT crystal with a cut angle of 40° 28 1 had very small change in resonance 
frequency when its temperature was changed by ±10°K while submerged in 
liquid helium. An oscillator circuit was also designed to operate at low 
temperatures in high vacuum. The oscillator was kept warmer than the walls 
by spot heaters installed in its chassis. The crystal was provided with two 
vacuum deposited gold electrodes which were connected to the oscillator by 
a coaxial cable long enough to absorb the temperature jump between the oscil- 
lator package and crystal. Its length, however, was limited by the noise 
introduced in the oscillator, since such a lead tends to act like an antenna. 
Based on this type of mounting, the frequency was stabilized in high vacuum 
and low temperatures to one part per 5 X 10 , which corresponds to a 
sensitivity of more than a molecular layer. 

In order to maintain the crystal temperatures as close as possible to 
the wall temperature, a U-shaped copper plate was adopted for the crystal 
chassis. Four crystals were mounted on the plate by two diagonally opposed, 
low outgassing, silver epoxy spots (see Figure 6a). The front side of the 
copper plate was provided with four holes of 1. 77 cm (0. 5 in. ) diameter to 
expose the area of the crystal that was covered by the electrode (see Fig- 
ure 6b). A complete set of four crystals with the block of remote oscillators 
can be seen in Figure 6c. The copper plate was provided with spot heaters 
and a germanium thermistor, allowing for the thermal control necessary for 
the testing. The thermal control will be described in the Procedures section. 
The plate was clamped around the gaseous helium tubing that refrigerates the 
molecular trap, using indium to improve the joint conductance. Two of these 
sets of crystals were prepared and installed in the locations indicated in 
Figure 7 and identified as crystal group I and II. Group III referred to in the 
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same figure corresponds to a different QCM design that will not be described 
in detail in this paper. It consists of a temperature compensated QCM that 
can be operated in a range of temperatures of more than 100°K. It is 
described in complete detail in Reference (8). Figure 7 also shows all the 
dimensions involved in the test set up and used in the reduction of the data 
for this paper. 


III. PROCEDURES 

The output of the crystals was fed into a frequency counter and recorded. 
The data systems utilized are the same as described in Reference (8) although 
no filter was necessary for this test series to clean the crystal output signals. 
The digitized signals were recorded on paper tape, and a different file was 
opened for each run. 

The variables considered in each run consisted of a nozzle type, its 
vertical position, kind of gas, and plenum pressure level. Prior to the run, 
the selected nozzle was warmed to room temperature by the coaxial heater 
and its temperature was stabilized by a slight flow of gas. The plenum pres- 
sure was then set to the desired level and the crystal data was recorded for 
about one to two minutes, at the end of which the gas valve was shut. Record- 
ing of the crystal outputs then continued for some time to check for possible 
desorption. Since the number of tests was so large, the sensitivity of the 
crystals was affected from time to time by the amount of mass deposited on 
them. It then became necessary to stop the test momentarily and warm up 
the crystals by energizing the heaters provided for each crystal package. In 
this manner a fast thermal desorption allowed for a quick and practical clean- 
up and recovery of the crystal, which was ready to operate again one hour 
after its heater was turned back to the off position. 

All the systems operated nominally and the total series of tests was 
accomplished in eight days. At no time was it necessary to pump out the 
chamber gases, although at the end of the test the background pressure had 
decayed to Z X 10 N/m 2 (1.5 X 10 torr) from an initial value of Z X 10“*® 
N/m 2 (1. 5 X 10“ 12 torr). 
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IV. RESULTS 


The results presented in this paper are summarized in Table II. Each 
of the five nozzles was tested with carbon dioxide and nitrogen, and each 
gas was run at several plenum pressures. Furthermore, each case was run 
for two or three different positions of the nozzle, varying the relative posi- 
tion of the plenum with respect to the QCMs. All nozzles were run at approxi- 
mately 294°K (70°F) plenum temperature. 

The cases described in Table II were analyzed by grouping the cases 
that differ only in the distance of the nozzle to the Molsink top. The raw data 
is presented as an array in which the frequency of each QCM is tabulated 
versus time. From this table and by using Equation 1, the following equation 
is obtained: 


M _ mas flux rate arriving at the crystal _ dlvl 
A c area of the crystal ~ dA c 

Furthermore, dA c can be expressed as 

dA c = r2 cos 4> 


where r is the distance from the nozzle exit section to the center of the crys- 
tal, <j> is the angle between the normal to the crystal and T, and df 2 is the 
solid angle substended by dA c from the center of the nozzle exit section. 
From this expression one obtains 

iCl r 2 

= mass flux per unit solid angle =5 7 — - 

c 6 A c cos4> 

where 6 is the angular distance from the crystal to the plume axis. If one 
normalizes this by its value at the axis, one obtains 




f(e) 
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The function f(0) has received considerable attention in the literature, and 
the data in this paper has been reduced to provide information about this 
function. The value of (dm/dD) 0 _ Q obtained by the Hill and Drapper scheme 
(Reference (9)) has been adopted, and the resulting data compared with the 
f( 0 ) proposed by those investigators, i. e. , 

f(0) = exp j-[\^(l “ cos 0)^lj 


where 


X 



C F /C FM AX )] 


-1 


is the plume slenderness coefficient, and 



is the mass normalization factor at the axis center line and m w - nozzle total 
mass flow rate. The data is being presented in Figures 8 to 17 in semi- 
logarithmic scale, i. e. , 


f(6) 



( 9 ) 


where is the mass coefficient. The normalization factor (dm/df2) 0 _ q 
is also given for convenience. Tables II to VII contain the numerical values 
of K m for several 0's. Because of the logarithmic characteristics of the 
ordinate, the cases corresponding to different pressures for the same nozzle- 
gas combination have been grouped together. 

V. DISCUSSION 

The results described in the above paragraphs constitute a representa- 
tive sampling of the total amounts of data obtained. In the majority of cases, 
the data is accurate within ±5% as measured by the QCMs. There are, how- 
ever, some cases in which the QCMs readings were relatively poor, and they 
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have been marked in Table II by an asterisk. Nevertheless, these cases 
were the exception and were mostly due to QCM saturation because of cumu- 
lative mass deposits and are easy to detect and identify. 

A question arises about the correlation between the mass flux as 
registered by the crystal and the actual mass flux impinging on the crystal 
surface. In order to check such a correlation, QCMs frequencies were 
compared before and after closing the gas valve. While the QCM presented 
a constantly decreasing frequency before closing the hand valve, the frequency 
was stabilized immediately to its last value once the hand valve was closed. 
This behavior indicates that all the mass arriving to the crystal was instantly 
captured by its surface, as was expected because of the low crystal tempera- 
ture as compared with the nitrogen condensation temperature. Migration 
effects, on the other hand, were responsible for dumping relatively large 
amounts of mass from the copper plate onto the crystals, but these rates 
were extremely slow as compared with the plume gas deposition, and they 
only caused saturation problems and signal modulations with periods on the 
order of hours. 

One peculiarity easily observed in Figures 8 and 9 is the fact that dif- 
ferent groupings of data are obtained for the same plume under different 
relative positions within the chamber and respect to QCMs. It is interesting 
to note that the intermediate grouping corresponds in both runs to the lowest 
position of the nozzle, i. e. , 60. 96 cm (Z. 0 ft. ) from the top, the upper group- 
ing corresponds to the intermediate position of the nozzle (45. 72 cm, 1.5 ft. ) 
and the lower one to the highest position (30. 28 cm, 1. 0 ft. ), The reason for 
this is not well understood as yet, but could be caused by different QCM sensi- 
tivities (quite improbable), by different residual recirculation effects (not 
consistent with the relative position of the groupings of data), and/or gas 
surface effects at the nozzle lip (these effects are randomly distributed and 
are inconsistent with the organized order encountered in Figures 8 and 9). 

This question will be left open at this time, although the data scatter does not 
appear too adverse. Because of the small expansion ratio of Nozzle 1, the 
location of the data approaches that region of the plume in which the boundary 
layer correction departure from the inviscid theory begins, which shows that 
the QCMs readings are consistent with the theory in the region of the 
theory 1 s validity. A comparison of Figures 10 and 14, and 11 and 15, 
shows the effects of a flat nozzle edge versus a sharp edge. It can be seen 
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that the nozzle with a sharp edge, as expected, spills more mass in the back 
region than the flat edged one. The same scattering mentioned before can be 
observed here. Note also that, in most cases, an increase in pressure 
reduces the back flow mass flux, although this effect seems to be less accen- 
tuated than the data spread due to the geometrical position of the plume in the 
chamber. 

In Figures 12 and 13 one can see that the high expansion ratio of the 
nozzle produces a very slender theoretical plume. However, the back flow 
mass flux is little affected by the high expansion rates and thus remains of 
the same order of magnitude as in the cases previously mentioned. Figures 
16 and 17 show the effect of the conical angle. Comparison with Figures 10 
and 11 does not indicate any outstanding difference between a 15° or a 25° 
cone angle, except in those regions closer to the axis where the data tends 
to increase or decrease with the cone angle. This variation is probably due 
to two dimensional effects at the exit plane. 

VI. SUMMARY 

Forty-four cases of back flow measurements have been reported. 
Although some data spread is observed, enough information is presented 
to show the order of magnitude involved in the correction of existing theories. 
This set of data is the most comprehensive ever obtained and provides an 
answer to the many expressed concerns with estimates of exhaust gases in 
those regions of the plume that depart more than 90° from the center line. 
This could be of significance if materials incompatible with the gaseous 
exhaust products are used in this vicinity, or if a science experiment con- 
ducted in these regions could be affected by the presence of these gases. 

VII. CONCLUSIONS 

From the previous discussion, one concludes: 

(1) The cryogenic QCMs appear to be the ideal type of instrument to 
perform plume far field mass flux measurements. 

(2) The back flow measurements seem to converge to the predicted 
theoretical values for smaller departures from the plume axis. 
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(3) 


The specific heat ratio and/or molecular weight of the gas does 
not seem to affect the data drastically. 

(4) A plenum pressure increase reduces the mass flux into the back 
flow region. 

(5) The flat edge acts as a molecular beam deflector; that is, it 
screens the back flow region from molecular bombardment. 

(6) The cone angle of the nozzle affects only that data which is rela- 
tively close to the axis. 

(7) A large expansion ratio does not affect the back flow mass flux. 

As a result of these conclusions, it is inferred that the mechanisms of 
gas-surface interaction that take place at the nozzle edge may be dictating 
the amount of mass flux expected for large departure angles (over 100°), 
while the boundary layer expansion may regulate the far field structure from 
40° to 100° 
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NOMENCLATURE 


C F 

C M 


C 


T 


f 

K M 

M 

M 


m(6) 

r 


r 

s 

\ 

<J> 

n 

e 


2 

= area of the crystal (cm ) 
= thrust coefficient 


= crystal mass coefficient 



= crystal temperature coefficient 


, (Hz cm 2 /g) 


- average crystal frequency, (MHz) 

= crystal frequency, (Hz) 

= back flow mass flux coefficient 
= mass deposit on crystal (g/cm^) 

= mass flux rate arriving at the crystal (g/cm^/sec) 

= mass flux rate at 0 location, (g/cm^/sec) 

= radius vector from the nozzle exit plane center to the 
crystal center 



(cm) 


= nozzle edge width (cm) 

= plume slenderness parameter 

= angle between the normal to the crystal and the radius 
vector r 

= solid angle 

= angular distance from plume axis 
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Table I. Geometric characteristics of the five tested nozzles 
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Table II. Summary of cases presented and analyzed 


Case No. 

Nozzle 

" 

Gas 

Pressure 
N/m^ (psia) 

Distance from Molsink 
top-cm (in. ) 

i 

i 

N 2 

2. 19 X 10 4 (3. 2) 

30. 48 (12) 

2 

i 

N 2 

2. 19 X 10 4 (3. 2) 

45. 72 (18) 

3 

i 

N 2 

2. 19 X 10 4 {3. 2) 

60. 96 (24) 

4 

i 

co 2 

2. 19 X 10 4 (3. 2) 

30. 48 (12) 

5* 

i 

co 2 

2. 19 X 10 4 (3. 2) 

45. 72 (18) 

6 

i 

co 2 

2. 19 X 10 4 (3.2) 

60. 96 (24) 

7* 

2 

N 2 

2. 19 X 10 4 (3.2) 

66. 04 (26) 

8* 

2 

N 2 

8. 5 X 10 4 (12. 4) 

45. 72 (18) 

9 

2 

N 2 

1. 2 X 10 5 (17. 5) 

45. 7 2 (18) 

10* 

2 

C °2 

2. 19 X 10 4 (3. 2) 

45. 72 (18) 

11 

2 

C °2 

2. 19 X 10 4 (3. 2) 

66. 04 (26) 

12 

2 

C °2 

8. 5 X 10 4 (12. 4) 

45. 72 (18) 

13* 

2 

C °2 

1. 2 X 10 5 (17. 5) 

45. 72 (18) 

14 

3 

N 2 

2. 19 X 10 4 (3.2) 

60. 96 (24) 

15 

3 

N 2 

8. 5 X 10 4 (12. 4) 

60. 96 (24) 

16 

3 

N 2 

1.2 X 10 5 (17. 5) 

60. 96 (24) 

17 

3 

co 2 

2. 19 X 10 4 (3. 2) 

60. 96 (24) 

18 

3 

C °2 

8. 5 X 10 4 (12. 4) 

60. 96 (24) 

19 

3 

C °2 

1. 7 X 10 5 (25. 0) 

60. 96 (24) 

20 

4 

N 2 

2. 19 X 10 4 (3. 2) 

45. 72 (18) 

21 

4 

N 2 

2. 19 X 10 4 (3. 2) 

66. 04 (26) 

22 

4 

N 2 

8. 5 X 10 4 (12. 4) 

45. 72 (18) 

23 

4 

N 2 

8. 5 X 10 4 (12. 4) 

66. 04 (26) 

*Cases suspected of containing some poor QCM readings 
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Table II (contd) 



Nozzle 

Gas 

Pressure 
N/m^ (psia) 

Distance from Molsink 
top-cm (in. ) 

4 

N 2 

8. 5 X 10 4 (12. 4) 

91. 44 (36) 

4 

N 2 

1. 2 X 10 5 (17. 5) 

45. 72 (18) 

4 

N 2 

1. 2 X 10 5 (17. 5) 

66. 04 (26) 

4 

N 2 

1. 2 X 10 5 (17. 5) 

91. 44 (36) 

4 

C °2 

2. 19 X 10 4 (3. 2) 

66. 04 (26) 

4 

co 2 

2. 19 X 10 4 (3. 2) 

91. 44 (36) 

4 

co 2 

8. 5 X 10 4 (12. 4) 

66. 04 (26) 

4 

co 2 

5 

1. 2 X 10 (17. 5) 

66. 04 (26) 

4 

co 2 

1. 2 X 10 5 (17. 5) 

91. 44 (36) 

5 

N 2 

2. 19 X 10 4 (3. 2) 

50.8 (20) 

5 

N 2 

2. 19 X 10 4 (3. 2) 

76.2 (30) 

5 

N 2 

8. 5 X 10 4 (12. 4) 

50.8 (20) 

5 

N 2 

8. 5 X 10 4 (12. 4) 

76.2 (30) 

5 

N 2 

1. 2 X 10 5 (17. 5) 

50.8 (20) 

5 

N 2 

1.2 X 10 5 (17. 5) 

76.2 (30) 

5, 

co 2 

2. 19 X 10 4 (3. 2) 

50.8 (20) 

5 

C °2 

2. 19 X 10 4 (3. 2) 

76.2 (30) 

5 

C °2 

8. 5 X 10 4 (12. 4) 

50.8 (20) 

5 

C °2 

8. 5 X 10 4 (12. 4) 

76.2 (30) 

5 

C °2 

1.2 X 10 5 (17. 5) 

50.8 (20) 

5 

C °2 

1. 2 X 10 5 (17. 5) 

76.2 (30) 


*Cases suspected of containing some poor QCM readings 
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Table III. Numerical values of back flow mass coefficient 
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Table III (contd) 


(S?)e /(^)e - o ' £W - 10 ' KM 

Carbon dioxide 

4 2 

Plenum pressure 2. 19 X 10 N/m (3.2 psia) 

(^)e = o =3 - 05g/sec 

0 

k m 

56. 97 

2. 20 

64. 33 

2. 72 

71. 43 

2. 83 

73. 21 

2. 92 

75. 09 

3. 00 

77. 06 

3. 09 

87. 34 

3. 51 

89. 65 

3. 57 

92. 16 

3. 67 

94. 89 

3. 76 

102. 32 

3. 73 

104. 72 

3. 80 

106. 57 

3. 31 

107. 16 

3. 85 

109. 42 

3. 37 

112. 40 

3. 46 

115. 52 

3. 54 

122. 42 

4. 27 

125. 42 

4. 31 


18 
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Table IV. Numerical values of back flow mass coefficient 

versus 0 for Nozzle 2 


M,M.o 

Nitrogen 

4 2 

Plenum pressure 2. 19 X 10 N/m (3. 2 psia) 

(te) d = 0 = °- 585 § /sec 

e 

K. , 


M 

67. 68 

2. 93 

103. 95 

4. 23 

106. 19 

4. 31 

108. 46 

4. 28 

110. 73 

4. 35 

125. 31 

4. 27 

127. 80 

4. 27 

132. 84 

4. 34 

(§)./(*?). . , ■ *> ■ 

Nitrogen 

4 2 

Plenum pressure 8. 5 X 10 N/m (12. 4 psia) 

(^)e = o = 2 ' 27g/sec 

0 

K. , 


M 

57. 81 

2. 71 

85. 63 

3. 29 

92. 15 

4. 00 

108. 22 

4. 54 

113. 88 

5. 46 
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Table IV (contd) 


/dm\ // dm\ 

We/W 


9 = 0 


= f( 0) = 10 


5 tvt 1 ^ 


Plenum pressure 1.2 X 10 N/m (17.5psia) 


113. 88 




= 10 


Nitrogen 




Plenum pressure 2. 19 X 10 N/m (3. 2 psia 


Carbon dioxide 


= 1.38 g/ sec 


57. 81 
67. 68 
85. 63 
92. 15 
103. 95 
105. 57 
106. 19 
108. 22 
108. 46 
110. 73 
113. 88 
125. 31 
127. 80 
132. 84 


2. 95 

3. 40 

4. 45 
4. 24 

4. 64 

5. 15 
4. 68 
4. 60 
4. 72 

4. 77 

5. 07 
4. 89 
4. 90 
4. 91 
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Table IV (contd) 


(t),/(t), . o ■ f(e > * 

Carbon dioxide 

4 2 

Plenum pressure 8.5 X 10 N/m (12. 4 psia) 

(t) 6 ,o =5 ' 34s/sec 

0 

k m 

57. 81 

o 

00 

rvi, 

92. 15 

4. 19 

108. 22 

4. 75 

110. 99 

5. 04 

(df) e /(df) 0 = o = f(e) = 10 M 

Carbon dioxide 

5 2 

Plenum pressure 1.2 X 10 N/m (17. 5 psia) 

(§) e = 0 

0 

k m 

57. 81 

2. 46 

110. 99 

5. 80 

00 

00 

ro 

6. 22 
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Table V. Numerical values of back flow mass coefficient 

K,, versus 0 for Nozzle 3 
M 


/dm\ // dm\ 

\ dQ /e/ \ dS2 /e = 


= f(0) = 10 


Nitrogen 


Plenum pressure 2. 19 X 10 N/m (3. 2 psia) 


= 0. 227 g/ sec 



/ drii\ //drrA 

\an)e/\d a) e = 


= f(e ) = 10 


4 , 2 


Nitrogen 


Plenum pressure 8.5 X 10 N/m (12.4 psia) 


(®)e- 


= 0. 879 g/sec 


66. 08 
118. 75 
121. 14 
126. 06 


3. 02 
4. 10 
4. 49 
4. 62 


/driA //dm\ - £(0) = 10 

We/\ dn /e = 0 


5 2 

Plenum pressure 1.2 XlO N/m (17. 5 psia) 


Nitrogen 


dm\ 

d^/e = 


= 1 . Z4 g/sec 
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Table V (contd) 
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mi mu 1 1 M iif HP I 1 I I II ill l 


Table VI. Numerical values of back flow mass coefficient 

versus 0 for Nozzle 4 


(i)e /(£)* . 0 ■ f,e > ■ io ' KM 

Nitrogen 

4 2 

|MH 


Plenum pressure 2. 19 X 10 N/m (3. 2 psia) 

■Enfl 

| = 0. 585 g/ sec 




9 

PL , 


M 

68. 63 

2. 84 

85. 93 

3. 53 

87. 89 

3. 60 

89. 91 

3. 67 

92. 00 

3. 76 

105. 14 

3. 67 

107. 25 

3. 70 

109. 18 

4. 34 

109. 38 

3. 74 

111. 80 

4. 31 

122. 89 

3. 83 

125. 21 

3. 85 

129. 95 

3. 88 



/ rTm\ //rim\ M 


(S?) /(sn) = f ( 0 ) = 10 

Nitrogen 

\ dn /e / \ dfi /e = o 


4 2 

Plenum pressure 8.5 X 10 N/m (12.4 psia) 

ms 

i = 2* 267 g/ sec 


\ u “/e = o 

0 

k m 

59. 03 

2. 89 

82. 57 

3. 42 

85. 93 

3. 60 

87. 89 

3. 67 

89. 91 

3. 74 

91. 99 

3. 80 

104. 27 

4. 29 

105. 13 

3.91 

106. 67 

4. 37 

107. 25 

3. 94 

109. 38 

3. 97 

121. 39 

4. 65 

122. 88 

4. 13 

123. 21 

4. 71 

125. 01 

4. 74 

125. 21 

4. 15 


JPL Technical Memorandum 33-620 


IT 


24 








Table VI (contd) 


= f( 9 ) = 10 


(jg)e/(^).-o ^ (9) = 1 ° 

Plenum pressure 8. 5 X 10 N/m (12. 4 psia) 


9 = 0 


126. 79 
129. 95 
138. 66 
140. 51 
144. 16 


/dm\ //dm\ = f(e) = 10 

\dn) e / \dn ) e = o 


4. 82 
4. 19 
4. 89 

4. 95 

5. 00 


Plenum pressure 1.2 X 10 N/m (17. 5 psia) _ 


68. 63 
82. 57 
85. 93 
87. 89 
89. 91 
92 . 00 

104. 27 

105. 13 

106. 67 
107. 25 

109. 38 

121. 39 
122. 89 

123. 21 
125. 01 

125. 21 
126. 69 
129. 95 

138. 67 
140. 51 
144. 16 


2 . 92 
3. 39 
3. 50 

3. 58 
3.65 
3.71 

4. 17 

3. 86 

4. 24 
3. 90 

3. 94 

4. 60 
4. 12 
4. 66 
4. 73 
4. 14 
4. 77 
4. 18 

4. 95 

5. 04 
5.09 
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Table VI (contd) 


(4/W..0 f(e) = 10""“ 

Carbon dioxide 

4 2 

Plenum pressure 2. 19 X 10 N /m (3. 2 psia) 

W.-o =i - 38g/sec 

0 

K 


M 

68. 63 

3. 18 

82. 57 

3. 72 

105. 14 

4. 04 

107. 25 

4. 06 

109. 38 

4. 09 

121. 39 

4. 98 

123. 21 

5. 10 

125. 01 


125. 21 

4. 25 

126. 79 

5. 25 

127. 57 

5. 31 

138. 67 

5. 35 

140. 51 

5. 59 

144. 17 

5. 57 

W. /(&)..„ = f(e) = io ' Km 

Carbon dioxide 

Plenum pressure 8. 5 X 10^ N/m^ (12. 4 psia) 

(if)e = o = 5 - 34 g/ sec 

0 

K r 


M 

68. 63 

3. 22 

105. 14 

4. 12 

107. 25 

4. 15 

109. 38 

4. 19 

122. 89 

4. 39 

125. 21 

L 

4. 41 
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Table VI (contd) 
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Table VII. Numerical values of back flow mass coefficient 

versus 0 for Nozzle 5 



Nitr ogen 

4 2 

Plenum pressure 2. 19 X 10 N/m (3. 2 psia) 

( TTTT' ) = 0.585 g/sec 

\ d ^/e = 0 B 

e 

K_ , 


M 

60. 11 

2. 68 

92. 46 

3. 92 

94. 70 

3. 98 

97. 00 

4. 09 

111.11 

4. 51 

112. 33 

3. 64 

113. 78 

4. 58 

114. 48 

3.71 

116. 63 

3. 72 

118. 77 

3. 76 

119. 42 

4. 72 

132. 71 

3.84 

134. 97 

3. 85 

139. 47 

3. 97 

(&)./(*).. o ^> = I0 ' KM 

Nitrogen 

4 2 

Plenum pressure 8. 5 X 10 N/m (12. 4 psia) 

fe) =2.27g/se c 


\ d ^/e = o 

0 

K, , 


M 

60. 11 

2. 95 

73. 24 

3. 22 

92. 46 

3. 96 

94. 70 

4. 03 

97. 00 

4. 13 

111.11 

4. 76 

113. 76 

4. 89 

114. 48 

4. 12 

116. 63 

4. 16 

118. 77 

4. 18 

119. 42 

5. 13 

132. 71 

4. 27 

134. 97 

4. 29 

139. 47 

4. 59 
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Table VII (contd) 


/dm\ //dm\ 

V 1 “/e/We 

- K M 

= f(0) = 10 

= 0 

Nitrogen 

Plenum pressure 

1. 2 X 10 N/m (17.5 psia) 

(£?)e = o = 3 ' 20g/sec 


e 

k m 


60. 11 

2. 88 


73. 24 

3. 27 


92. 46 

3. 89 


94. 70 

3. 97 


97. 00 

4. 05 


111.11 

4. 69 


113. 78 

4. 81 


114. 48 

4. 20 


116. 63 

4. 23 


118. 77 

4. 25 


119. 42 

5. 08 


132. 71 

4. 37 


134. 97 

4. 38 


139. 47 

4. 88 

/ dm\ //dm\ 

\ d ^ /e /\dfi ye 

-K 

= f(0) = 10 

= 0 

Carbon dioxide 

Plenum pressure 

4 2 

2. 19 X 10 N/m (3. 2 psia) 

= 1.38 g/sec 

\d«y 0 = 0 B 


0 

k m 


60. 11 

3. 14 


73. 24 

3. 28 


92. 46 

4. 19 


94. 70 

4. 29 


97.00 

4. 38 


111.11 

4. 93 


112. 33 

4. 16 


113. 78 

5. 05 


114. 48 

4. 18 


116. 63 

4. 21 


118. 77 

4. 23 


119. 42 

5. 22 


132. 71 

4. 41 


134. 97 

4. 34 


137. 22 

4. 43 


139. 47 

4. 52 
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Table VII (contd) 


fdrn\ //dm \ 

/e/VWe = 0 


= f( 0 ) = 10 


Carbon dioxide 


Plenum pressure 


8. 5 X 10^ N/m^ (12. 4 psia) 



= 5.34 g/ sec 


60. 11 
73. 24 
92. 46 
94. 70 
97. 00 
111.11 
112. 33 
113. 78 
114. 48 
116. 63 
118. 77 
119. 42 
134. 97 
139. 47 


3. 17 
3. 40 

4. 22 
4. 30 
4. 38 
4. 99 

4. 40 

5. 11 
4. 44 
4. 49 

4. 51 

5. 33 
4. 66 
4. 69 



Carbon dioxide 


= 7. 53 g/ sec 


2 . 8 ‘ 
3. 4( 


4. 00 
4. 09 
4. 16 
4. 72 
4. 35 
4. 83 
4. 40 
4. 44 

4. 47 

5. 04 
4. 60 
4. 69 
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Fig. 1. Molsink vacuum chamber and auxiliary equipment 






Arrangement of the 5 plenum-nozzle assembly units in the Molsink 
(all dimensions are in cm) 
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Af = C M AM + C T AT 

2 

AM (gm/cm ) 

C M = <*>( T, CUT ANGLES) 

Af (Hz) 

C T = *(T, CUT ANGLES) 

F (MHz) 
c 

If C, = 0, AM = 4.35 x l(f 7 

1 

c 




Fig. 5. Quartz crystal, electrodes, and demonstration of the thickness 
vibration mode 
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CHASSIS HEATER 



36 


Reari;view of 4 cryogenic crystals and other components mounted on a common chassis 
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Fig. 6c. Crystal chassis and the remote oscillators block 
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CRYSTAL GROUP A 


n 

]_!" 7 


iMi 

>4 

u 


t 

2.54 

2.54 

2.54 

17.7 

2.54 



1.270 DIA HOLES 
W/CRYSTALS 
CENTERED A 

t 


4.12 RAD - 5 HOLES 
SPACED EQUALLY 


TOP VIEW 


NOTE: 

ALL QCM's AND 
NOZZLE No. 1 ARE 
CENTERED ON 
SECTION PLANE 


CRYSTAL GROUP B 


DETAIL A 


ALL DIMENSIONS TAKEN 
FROM PLUMB LINE TO 
CENTER OF CRYSTAL 
GROUPINGS 

1.270 DIA HOLES 
W/CRYSTALS 
CENTERED 


^ 1.42 


1 ( 
5.08 

i ' 


V 

* 

rh 

i ■ 


T 

2.2 

!2 

k3.81* 

-•—5.08 


DETAIL B 


♦ASTERISK IS LOCATION OF CENTER OF 
CRYSTAL GROUPING - DIMENSIONS 
LOCATING CRYSTAL GROUPS IN 
MOLSINK ARE TAKEN FROM ASTERISKS 
TO PLUMB LINE 



Fig. 7. Location and position of the QCM units in the Molsink chamber 
(all dimensions are in cm) 
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MASS FLUX PER UNIT SOLID ANGLE 



Fig. 8. Mass flux measurement data reduced and compared with the Hill and 
Draper approximation for Nozzle 1 and N 2 gas 
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MASS FLUX PER UNIT SOLID ANGLE 



Fig. 9. Mass flux measurement data reduced and compared with the Hill and 
Draper approximation for Nozzle 1 and CCX> gas 
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ANGULAR DISTANCE FROM AXIS (0) deg 

Fig. 10, Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 2 and N£ gas 
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FLUX PER UNIT SOLID ANGLE (f£)/(S 



ANGULAR DISTANCE FROM AXIS (. 9 ), deg 


Fig. 12. Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 3 and gas 
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MASS FLUX PER UNIT SOLID ANGLE 


PLENUM 



PRESSURE 

(dm/dft) e _ 0 


N/m^ 

(PSIA) 


o 

2.19 x 10 4 

(3.2) 

0.61 

X 

6.8 x 10 4 

(10.0) 

2.36 

n 

1.7 x 10 5 

(2 5) 

3.33 


25° 


1 — 0.127 cm (0.05 in.) 
240 


NOZZLE 

GAS 

TEMP 


V 2 

294 K 


ANGULAR DISTANCE FROM AXIS (60, deg 

Fig. 13. Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 3 and CO^ gas 
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MASS FLUX PER UNIT SOLIO ANGLE 



ANGULAR DISTANCE FROM AXIS W), deg 

Fig. 14. Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 4 and gas 
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MASS FLUX PER UNIT SOLID ANGLE 



Fig. 15. Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 4 and CO^ gas 
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MASS FLUX PER UNIT SOLID ANGLE 


PLENUM 



PRESSURE 



N/m 2 

(PSIA) 


o 

2.19 x 10 4 

(3.2) 

0.585 

X 

8.5 x 10 4 

(12.4) 

2.267 

□ 

1.2 x 10 5 

(17.5) 

3.199 


15° 


>— 0.254 cm (0,1 in. 
6 = 60 


v-/ 

%° 


NOZZLE 

GAS 

TEMP 


5 

N2 

294 K 


ANGULAR DISTANCE FROM AXIS (0), deg 

Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 5 and gas 
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MASS FLUX PER UNIT SOLID ANGLE 


PLENUM 

PRESSURE 


(dm/dfl^Q ci/soc 



N/m^ 

CPSIA) 


o 

2.19 x 10 4 

(3.2) 

1.37 

X 

8.5 x 10 4 

(12.4) 

5.33 

□ 

1.2 x 10 5 

(17.5) 

7.53 


15° \ d 


1 — 0.254 cm (0. 1 in. ) 
€ = 60 


NOZZLE 

GAS 

TEMP 


^^2 
294 K 


ANGULAR DISTANCE FROM AXIS (9), deg 

Fig. 17. Mass flux measurement data reduced and compared with the Hill 
and Draper approximation for Nozzle 5 and CO^ gas 
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